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Abstract 

Simulations of the Greenland Ice Sheet are carried out with a high-resolution 
version of the ice-sheet model SICOPOLIS for several global- warming scenarios for 
the period 1990-2350. In particular, the impact of surface-meltwater-induced accel- 
eration of basal sliding on the stability of the ice sheet is investigated. A parame- 
terization for the acceleration effect is developed for which modelled and measured 
mass losses of the ice sheet in the early 21st century agree well. The main findings 
of the simulations are: (i) the ice sheet is generally very susceptible to global warm- 
ing on time-scales of centuries, (ii) surface-meltwater-induced acceleration of basal 
sliding leads to a pronounced speed-up of ice streams and outlet glaciers, and (iii) 
this ice-dynamical effect accelerates the decay of the Greenland Ice Sheet as a whole 
significantly, but not catastrophically, in the 21st century and beyond. 



1 Introduction 



In Chapter 10 ("Global Climate Projections") of the Fourth Assessment Report (AR4) of 
the United Nations Intergovernmental Panel on Climate Change (IPCC), an increase of 
the mean global sea level by 18-59 cm for the 21st century (more precisely: 2090-2099 
relative to 1980-1999) is projected for the six SRES marker scenarios Bl, B2, AlB, AIT, 



A2 and AlFI (Meehl et al. 2007). The main causes for this sea level rise are thermal 



expansion of sea water and melting of glaciers and small ice caps, and to a lesser extent 
changes of the surface mass balance of the Greenland and Antarctic Ice Sheets. However, 
recent observations suggest that ice flow dynamics could lead to additional sea level rise, 
and this problem is explicitly stated in the AR4: 

"Dynamical processes related to ice flow not included in current models hut 
suggested by recent observations could increase the vulnerability of the ice sheets 
to warming, increasing future sea level rise. Understanding of these processes 



is limited and there is no consensus on their magnitude." (IPCC 2007) 
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These conjectured dynamical processes are (i) basal sliding accelerated by surface meltwa- 
ter, (ii) reduced buttressing due to the loss of ice shelves, and (iii) penetration of ocean 
water under the ice. The first process is more relevant for the Greenland Ice Sheet (the 
focus of this study), whereas the latter two may affect the stability of the West Antarc- 
tic Ice Sheet. On the observational side, recent results from satellite gravity measure- 
ments for the period 2002-2005 indicate surprisingly large mass losses of 239 ± 23 km^ a~^ 



(0.66 lb 0.06 mm a ^ sea level equivalent) for the Greenland Ice Sheet (Chen et al. 2006). 



Furthermore, major outlet glaciers (Jacobshavn ice stream, Kangerdlugssuaq and Helheim 



glaciers) have sped up drastically since the 1990's (Rignot and Kanagaratnam 2006). 



2 Ice-sheet model SICOPOLIS 



For this study, we use the ice-sheet model SICOPOLIS ("Simulation COde for POLy- 
thermal Ice Sheets"), which simulates the large-scale dynamics and thermodynamics (ice 
extent, thickness, velocity, temperature, water content and age) of ice sheets three-dimen- 
sionally and as a function of time ( Greve|[l997 http://sicopolis.greveweb.net/). It is 
based on the shallow-ice approximation (e.g. Hutter 1983) and the rheology of an incom- 



pressible, heat-conducting power-law fluid [Glen's flow law, see Paterson (1994)]. External 
forcing is specified by (i) the mean annual air temperature at the ice surface, (ii) the surface 
mass balance (precipitation minus runoff), (iii) the sea level surrounding the ice sheet and 
(iv) the geothermal heat flux prescribed at the bottom of a lithospheric thermal bound- 
ary layer of 5 km thickness. For all simulations of this study, the horizontal resolution is 
10 km, the vertical resolution is 81 grid points for the cold-ice column, 11 grid points for 
the basal layer of temperate ice (if existing) and 11 grid points for the lithosphere layer, 
the time-step is 0.25 a, and the geothermal heat-flux distribution and model parameters 



are those used by Greve (2005). 



3 WREIOOO scenario 



Future global warming shall be prescribed exemplarily by the WREIOOO scenario, which 
assumes stabilisation of the atmospheric CO2 concentration at 1000 ppm (Cubasch et al. 



2001). The corresponding temperature change from 1990 (the "present") until 2350 is 
shown in Fig. [l] (along with similar scenarios with lower stabilisation concentrations). 

In order to obtain temperature and precipitation forcings for the Greenland Ice Sheet, 
the argumentation by Greve (2004) is followed. The surface temperatures shown in Fig. [l] 
are amplified by a factor 2 and imposed as uniform increases over the ice sheet, and the 
precipitations are assumed to increase by 5% per degree of ice-sheet-surface-temperature 
change. Surface melting is parameterized by the degree-day method in the version by 



Greve (2005). This approach is a critical simplification, and it should rather be replaced 



by an energy-balance model for more accurate results. However, since the objective of this 
study is to assess the impact of ice-dynamical processes on the decay of the Greenland Ice 
Sheet rather than making precise predictions of the decay itself, the use of the degree-day 
method is a reasonable compromise. 
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Global mean temperature change ATg for the scenarios WRE450, WRE550, 
WRE750 and WREIOOO (stabilisation scenarios for atmospheric CO2), by 



Figure 1: 
WRE650, 

Cubasch et al. (|2001|). In this study, only the WREIOOO scenario will be considered. 



4 Basal sliding 



Basal sliding is described by a Weertman-type sliding law in the form of Greve and Otsu 



and Le Meur 2001), 



(2007), based on Greve et al. (1998) and modified to allow for sub-melt sliding (Hindmarsh 



Vh = -Cb e 



X 



where fb is the basal-sliding velocity, Cb the sliding coefficient, Tb the basal shear traction in 
the bed plane, p the ice density, g the gravity acceleration, H the ice thickness and P\^=pgH 
the overburden pressure. The term e^b'''^™^ represents the exponentially diminishing sub- 
melt sliding, where is the temperature relative to pressure melting (in °C) and 7sms=l°C 
the sub-melt-sliding coefficient. 

Acceleration of basal sliding by surface meltwater is parameterized by an extension of 



the approach by Greve and Otsu (2007). The sliding coefficient is expressed as 



<-^b — <-^b 



1 + T^M' 



(2) 



where 



11.2 ma^^ Pa^^, M is the surface melt rate (runoff), 7 is the surface meltwater 
coefficient, and r and s are adjustable exponents. The idea behind this parameterization is 
to relate the sliding speed-up to the local surface melt rate, and account for the less efficient 
percolation of meltwater to the base in regions where the ice is thick by the dependency 
on the inverse ice thickness. 

The parameterization employed by Greve and Otsu (2007) corresponds to (r, s) = 



(0, 1). For this case, the authors show in their Appendix A that data reported by Zwally 



et al. ( |2002 ) from the Swiss Camp in central west Greenland give rise to the estimate 
7 = 0.1 a m~^. We will also consider the cases (r, s) = (0, 2) and (1, 1), for which the same 
arguments lead to estimates of 7 = 0.05 a^ m~^ and 100 a, respectively. 
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5 Simulations 



5.1 Set-up 



Five simulations with different settings for the acceleration of basal sliding by surface melt- 
water will be discussed in order to investigate to what extent this process can increase the 
vulnerability of the Greenland Ice Sheet to future warming. In run #1, acceleration of 
basal sliding by surface meltwater is not considered (7 = 0). Runs #2-4 have been con- 
ducted with (r, s) = (0, 1), (0, 2) and (1, 1), respectively, and values of 7 chosen according 
to the estimates given at the end of Sect. |4] (designated in Table [T| as "1 00%"). Run 7^5 
corresponds to the most extreme scenario considered by Greve and Otsu (2007), with the 
settings (r, s) = (0, 1) and 7 = 5 a m^^ (50 times the above estimate, therefore designated 
in Table [1] as "5000%"). All simulations start with the present-day ice sheet as initial 
condition, and the model time is from 1990 until 2350. 



Run 


7 


(r,s) 


V2002-2005 














[km^ a"-*-] 


[m SLE] 


[m SLE] 


[m SLE] 


#1 







37.7 


0.12 


0.55 


1.21 


#2 


100% 


(0,1) 


111.6 


0.14 


0.60 


1.31 


#3 


100% 


(0,2) 


172.0 


0.17 


0.68 


1.48 


#4 


100% 


(1,1) 


248.9 


0.18 


0.67 


1.42 


#5 


5000% 


(0,1) 


1627.8 


0.58 


1.51 


2.71 



Table 1: Set-up and results of runs #1-5. For the meaning of the basal-sliding parameters 
7, r, s see Eq. ^ and the accompanying text. V2002-2005 denotes the average loss of ice 
volume between 2002 and 2005, whereas AV2100, AV2200 and AV2300 are the losses of ice 
volume by 2100, 2200 and 2300, respectively, compared to 1990. The latter are expressed 
in meters of sea- level equivalent. 



5.2 Results 

An overview of the main results is given in Table [TJ The average loss of ice volume be- 



tween 2002 and 2005 can be compared with the measured value by Chen et al. (2006) of 
239 ± 23 km^a~^ (see introduction), thus providing an observational constraint for the 
simulations. Evidently, the ice- volume loss is far too small for run #1 (no acceleration 
of basal sliding by surface meltwater) , and it is too small by about a factor 2 for run #2 
[(r, s) = (0,1)]. By contrast, the agreement is quite good for run #3 [(r, s) = (0,2)] and 
very good for run #4 [(r, s) = (1,1)]. On the other hand, the extreme case of run #5 
[(r, s) = (0,1), very large 7] produces more than 6 times more ice- volume loss than ob- 
served. Therefore, runs #3 and 4 seem to be most realistic. 

From a theoretical point of view, the set-up of run #4 is preferable to that of run #3, 
because it is clear that the percolation of surface meltwater to the base will be the less 
efficient the thicker the ice is. This is accounted for in run #4, for which the acceleration 
of basal sliding decreases with increasing ice thickness (r = 1), whereas this is not the case 
in run #3 (r = 0). Consequently, run #4 shall be considered as the "best" simulation. 
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Comparison of the results of run #4 and run 7^1 (no acceleration of basal sliding by 
surface meltwater) shows that the contribution to sea-level rise by 2100 is ~ 50% larger 
for run #4 (0.18 vs. 0.12 m). The impact of the acceleration effect on ice flow becomes 
evident by inspection of Fig. |2] which shows the simulated surface velocities in 2100 for the 
two runs. Therefore, the acceleration of basal sliding by surface meltwater, which is most 
likely the major ice-dynamical process relevant for the Greenland Ice Sheet in the context 
of global warming, has a significant, but not catastrophic effect on the decay of the ice 
sheet in the 21st century. 

The absolute difference between the two runs becomes larger in the more distant future; 
however, the relative difference becomes smaller: by 2200 the contribution to sea-level rise 
is 0.12 m (~ 22%) larger for run #4, and by 2300 it is 0.21 m (~ 17%) larger. Figure g 
shows the simulated surface topographies in 2350 (at the end of the simulations). It is 
nicely illustrated that for both runs #1 and #4 the ice sheet shows a strong response on 
the imposed warming scenario and retreats all around the margin (most pronounced in the 
south-west), while the surf ace- meltwater- induced acceleration of basal sliding accounted 
for in run #4 speeds up the decay. 

Two additional simulations with larger exponents, namely (r, s) = (1, 2) and (2, 1), and 
values of 7 chosen in analogy to the "100%" runs #2-4, have also been conducted. For 
these cases, maximum surface velocities of more than 100 kma^^ occur close to the ice 
margin, which is unrealistic. Apparently, the speed-up effect is too pronounced for these 
settings, and so they have been discarded. 
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2100, run #1 2100, run #4 




X (km) X (km) 

Figure 2: Simulated surface velocity of the Greenland Ice Sheet. Left: Run #1 (no 
acceleration of basal sliding by surface meltwater), year 2100. Right: Run 7^4 ("best" 
simulation, including acceleration of basal sliding by surface meltwater), year 2100. Both 
simulations reproduce the organization of the drainage pattern into ice streams and outlet 
glaciers, despite the use of the shallow-ice approximation and the large-scale (10 km) 
resolution. The dynamic acceleration effect in run #4 is clearly visible all around the ice 
margin and leads to faster decay. 
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Figure 3: Simulated surface topography of the Greenland Ice Sheet. Top: Initial condition 
for the year 1990. Bottom left: Run ^1 (no acceleration of basal sliding by surface 
meltwater), year 2350. Bottom right: Run ^4 ("best" simulation, including acceleration 
of basal sliding by surface meltwater), year 2350. The difference between 1990 and 2350 
amounts to 1.59 m SLE for run ^1 and to 1.84 m SLE for run #4. 
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6 Conclusion 



The simulations discussed in this study suggest that ice-dynamical processes can speed up 
the decay of the Greenland Ice Sheet significantly in the 21st century and beyond. How- 
ever, a catastrophically accelerated decay can only be obtained with unrealistic parameter 
settings and thus seems to be unlikely. 
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